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Abstract 

The Horodecki family employed the Jaynes maximum-entropy principle, fitting the mean (bi) 

O : 

of the Bell-CHSH observable (B). This model was extended by Rajagopal by incorporating the 



dispersion (erf) of the observable, and by Canosa and Rossignoli, by generalizing the observable 
(B a ). We further extend the Horodecki one-parameter model in both these manners, obtaining a 



> 
CO 

o 

f- — ' i/iree-parameter (bi,o\,ct) two-qubit model, for which we find a highly interesting/intricate con- 

o ' 

I/-) . tinuum (— oo < a < oo) of Hilbert-Schmidt (HS) separability probabilities — in which, the golden 

o . 

ratio is featured. Our model can be contrasted with the three-parameter (b q , o~q, q) one of Abe and 
Rajagopal, which employs a (7(Tsallis)-parameter rather than a, and has simply q-invariant HS 
separability probabilities of |. Our results emerge in a study initially focused on embedding cer- 
tain information metrics over the two-level quantum systems into a q-framework. We find evidence 
that Srednicki's recently-stated biasedness criterion for noninformative priors yields rankings of 

X' 

5_i ■ priors fully consistent with an information-theoretic test of Clarke, previously applied to quantum 

systems by Slater. 
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I. INTRODUCTION 

Both Rajagopal [1], as well as Canosa and Rossignoli Q] have extended a well-known 
maximum-entropy model of the Horodecki family |3] to two-parameter models, but in dif- 
ferent fashions. Rajagopal incorporated the dispersion (o~f) of the Bell-CHSH observable 
(B) ^y], the mean (bi) of which is already fitted in the Horodecki model, while Canosa and 
Rossignoli fitted the mean of generalized Bell-CHSH observables (B a ). We combine their two 
approaches into a t/iree-parameter (61, af, a) model, for which we uncover a very interesting 
continuum (—00 < a < 00) of exact Hilbert-Schmidt separability probabilities (sec. IVH B 21 
Fig- HU (EH)) — in which, among other features, the golden ratio P, (| appears. 

Our model can be interestingly contrasted with a three-parameter (b q ,ag,q) one also 
of Abe and Rajagopal 7] (sec. IVH B|) . which incorporates the g-parameter (nonextensi- 
tity/Tsallis index/escort parameter), rather than the a-parameter of B a The contin- 

uum (over q) of separability probabilities {independently of the metric employed) is simply a 
constant (equal to | in the Hilbert-Schmidt case, and to the "silver mean" , y/2 - 1, for the 
Bures and other monotone metrics). We examine a two-parameter (61, a) Canosa- Rossignoli- 
type model also exhibiting an interesting (non-flat) continuum fsec. IV11 B 31 Fig. O (|35jl ). 
Additionally, we obtain exact (Hilbert-Schmidt and Bures) separability probabilities for the 
three-parameter Tsallis-Lloyd-Baranger model [8] (sec. IVII C|) 

Our results emerge in a study initially focused on embedding certain information metrics 
over the two-level quantum systems into a g-framework. We find evidence (sec. Fig. ^) 
that Srednicki's recently-stated biasedness criterion for noninformative priors ]9;j yields rank- 
ings of priors fully consistent with an information-theoretic test of Clarke [10], previously 
applied to quantum systems by Slater 

II. NONINFORMATIVITY OF PRIORS 

Some fifteen years ago, Wootters asserted that "there does not seem to be any natural 
measure on the set of all mixed states" [12I p. 1375]. He did, however, consider random 
density matrices with all eigenvalues fixed. He remarked that once "the eigenvalues are fixed, 
then all the density matrices in the ensemble are related to each other by the unitary group, 
so it is natural to use the unique unitarily invariant measure to define the ensemble" 



7 group, 
" Q, P. 
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1375] (cf. Q). 

Arguing somewhat similarly, Srednicki recently proposed that in choosing a prior distri- 
bution over density matrices, "we can use the principle of indifference, applied to the unitary 
symmetry of Hilbert space, to reduce the problem to one of choosing a probability distribu- 
tion for the eigenvalues of p. There is, however, no compelling rationale for any particular 
choice; in particular, we must decide how biased we are towards pure states" (3, p. 6]. 

To be specific, we find, in an analysis involving four prior probabilities (p's), that the 
information-theoretic-based comparative noninformativity test devised by Clarke yields a 
ranking 

PF q=1 > PB > PB q=1 trunc > PF (1) 

fully consistent (Fig. [TJ) with Srednicki's recently-stated criterion for priors of "biasedness to 
pure states" . Two of the priors are formed by extending certain metrics of quantum-theoretic 
interest from three- to four-dimensions — by incorporating the g-parameter. The three- 
dimensional metrics are the Bures (minimal monotone) metric over the two-level quantum 
systems and the Fisher information metric over the corresponding family of Husimi distribu- 
tions. The priors pb and pp are the (independent-of-q) normalized volume elements of these 
metrics, and Pf q=1 is the normalized volume element of the q- extended Fisher information 
metric, with q set to 1. While originally intended to similarly be the g-extension of the 
Bures metric, with q then set to 1, the prior PB q=1 trunc, actually entails the truncation of the 
only o/j^diagonal entry of the extended Bures metric tensor. Without this truncation, the 
(/-extended Bures volume element is null, as is also the case in three /wg/ier-dimensional quan- 
tum scenarios we examine (including the Abe-Rajagopal two-qubit states (sec. IV11 BJ) . for 
which we further find q-invariant Bures separability probabilities equal to — 1 [the "silver 
mean"] and Hilbert-Schmidt ones equal to | (sec. IVHB 1|) . and the Tsallis-Lloyd-Baranger 
two-qubit states (sec. IVIlT^l ). 



III. BURES METRIC 



The Bures (minimal monotone) metric — the volume element of which we normalize to 
obtain one (pb) of the four prior probability distributions of principal interest here — yields 
the statistical distance between neighboring mixed quantum states (p) 



14. 



15| | . It provides 



an embedding of the Fubini-Study metric \\Ji sec. IV], which gives the statistical distance 
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FIG. 1: Four univariate marginal prior probability distributions in the near-to-pure-state region 
r G [1 — e, 1], where r is the radial coordinate in the Bloch sphere representation of two-level 
quantum systems, and r = 1 corresponds to a pure state. The order of dominance fully complies 
with that (^Q) obtained by the information-theoretic-based comparative noninformativity test 

between neighboring pure quantum states (\ip)) (cf. ^3)- Hiibner gave an explicit formula 
for the Bures distance Q, p. 240] (cf. Q), 

d B (p 1 ,p 2 ) 2 = 2-2tr(p 1 1 /2 p 2 p 1 1 / Y 2 - (2) 



(3) 



Further, he expressed it in infinitesimal form as [18. eq. (10)] 



d B (p,p + dp) 2 = Si 



1 1 < i\dp\j > 



2 \i + Aj 

where the Aj's are the eigenvalues and the |i)'s, the eigenvectors of p. 



A. Three-Dimensional Case 



In |20j , using the familiar Bloch sphere (unit ball in Euclidean 3-space) representation of 
the two- level quantum systems (2x2 density matrices), 

P=U 1 + ZX ~ iy \, r 2 = x 2 + y 2 + z 2 <l, (4) 
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x + iy 1 — z 



it was found (cf. p. 128]), here converting from cartesian to spherical coordinates, 

x — r cos 9i, y — r sin 0\ cos 9 2 , z = r sin 6i sin 9 2 , (5) 

that 

rlr,(n n4-rln\ 2 = -( 

.(1-r 2 ) 



d B (p,p + dp) 2 = ^( q ~ r2 x dr 2 + dn 2 ) ■ (6) 



The term dr 2 corresponds to the radial component of the metric and dn 2 , the tangential 
component (dn 2 = r 2 dQ\ + r 2 sin 2 62). In the setting of the quantum monotone metrics - 
the Bures metric serving as the minimal monotone one — it is appropriate to express the 



tangential component of the Bures metric (JUj) in the form lg, eq. (3.17)], 



1 — r \ ~! 

: i + r )/B( TT ^)) , (7) 



where = ^ is an operator monotone function 21 1. 



The volume element of the Bures metric (jZJ is |?fz^j , which can be normalized to a prior 
probability distribution over the Bloch sphere, 

r 2 sin#i 
71^(1 — r z ) 

B. Four-Dimensional Case 

Now, we can construct a /bwr-dimensional family of (properly normalized/unit trace) 2x2 
escort density matrices (cf. [22]), 

/ \-i I 1 + z x — iy \ 

P {q} =((l-r) q + (l + ry) , (9) 

\x + iyl — zl 

for which q — 1 recovers the standard Bloch sphere representation (j3J). Applying Hiibner's 
formula ([3*]). we have found that the extended Bures metric (now incorporating the q- 
parameter) has the form 

<W„ (p,p + dp) 2 = + * (W< log 2 Wdq 2 + 4i ^ 1 ° g 1 W dqdr+ (10) 
+4 * 2W ° Ar* + { - 1+W " ? dn> 

^ ( r 2 _ X)2 ^ r 2 

where IV = that is, the ratio of the two eigenvalues of p. 

The tangential component of the metric (jlO)) can be expressed as ((1 + r) j Bures q (W)) , 
where 

f (t) - 2 ( 1 + ^ 1 + tg ) 2 (n) 

JBures q yt) — _|_ £g\2 ' ^ > 

This bivariate function appears (Fig. |2J) to be monotonically-increasing for any fixed q (cf. 

Q). 
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FIG. 2: The function fBures q (t) that yields the tangential component IjllJI of the extended (four- 
dimensional) Bures metric ()10|) 

Now, in the earlier stage of our analyses, due to a programming oversight, we were under 
the impression that the off-diagonal dqdr term of (fTUj) was simply zero. If we do employ 
the fully correct form, with this dqdr term included, we find that the volume element is 
null. This, of course, could not yield a meaningful prior probability distribution. However, 
having proceeded under the impression that the dqdr term was null, we obtained a number 
of results that appear to be of interest and of some relevance. Therefore, for much of this 
study, we will treat the dqdr term as null, and thus deal with a truncated g-extended Bures 
metric. 

In the context of the harmonic oscillators states, Pennini and Plastino have argued that, 
in addition to the physical lower bound (ignorance-amount) of q > that in a quantal 

n n 

regime, q can be no less than 1 |2j| - due to the Lieb bound on the Wehrl entropy |24|. 



However, for the two-level quantum systems to t 



the lower bound on the Wehrl entropy is \ 25t eq. (12)]. We, thus, consider q G [|,oo 



le study of which we restrict ourselves here, 



to be the range of possible values of the escort parameter q. In practice, though, we will, 
for numerical and graphical purposes and normalization of the (divergent over q G [1/2, oo]) 
truncated extended Bures volume element (Sec. IVIjl . consider that q G [|,500]. 

In Fig. El we show the two-dimensional marginal volume element of (fTUj) (after omission 
of the dqdr term) - - integrating out the spherical angles, 61,82, and leaving the radial 
coordinate r and the escort parameter q. In Fig. HJ further integrating out r, we show the 
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FIG. 3: Two-dimensional marginal of the truncated four-dimensional extended Bures volume ele- 
ment ([T0|) 
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FIG. 4: One-dimensional marginal ()12|) over q of the four-dimensional truncated extended Bures 
volume element (|10jl 

corresponding one-dimensional marginal volume element of (after omission of the dqdr 
term) over q. This (Fig. 0} has the exact expression 

tt(1 + log4) 



24g 



(12) 



This prior, thus, conforms to Jeffreys' rule — as opposed to the Bayes-Laplace rule, which 
would give a constant prior . 

In Fig. we integrate out q G [^,500], leaving a (deep bowl-shaped) one-dimensional 
marginal over r G [0,1]. (The corresponding marginal in the unextended Bures case is 
2(i_ r 2) ; so it is simply increasing with r, in that case.) The associated indefinite integral 
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FIG. 5: One-dimensional marginal (obtained from (|13[)) over r of the four-dimensional extended 
Bures volume element (|10|) after omission of the off-diagonal dqdr term 



over q is 



7r (q W q (3 + W 2q ) \og{W) - (1 + W q ) (2 + (1 + W q f log(l + W™))) 



(13) 



6 (-1+r 2 ) (1 + W q f log(W) 

(So, we obtain the function plotted in Fig. El by substituting g = 500 and q = ~ into (fT3*|) 
and taking the difference.) 

For g = 1, the extended Bures metric ffTTT^) reduces to 

o?ss U r es , =1 (p, P + dp) 2 = — (1 - r 2 ) log 2 Wdg 2 - - log Wdgdr + ds B (p, p + dp) 2 . (14) 

Normalizing the volume element of this metric — but first nullifying the off-diagonal dqdr 
term — to a (non-null) prior probability distribution over the Bloch sphere, we obtain (cf. 

©), 

3 r 2 sinfli log^ 

PB q=1 trunc- 4 ^+^4) > ^ 

one of the four priors that we rank (Fig.QJand (0)) both by the comparative noninformativity 
test and Srednicki's biasedness criterion. 



C. Comparative Noninformativities in the Bures Setting 



The relative entropy (Kullback-Leibler distance /information gain 

123, 28]) of pb with 



respect to PB q=1 trunc [which we denote Skl{pb-,Pb =1 trunc)] — that is, the expected value with 

respect to ps of log — is 0.101846 "nats" of information. Now, reversing arguments, 

SKL(j>B q= itrunc, Pb) = 0.0661775. (We use the natural logarithm, and not 2 as a base, with 
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one nat equalling 0.531 bits.) Let us convert - - using Bayes' rule — these two (prior) 
probability distributions to posterior probability distributions (posts and post BureSq=1 ), by 
assuming three pairs of spin measurements, one each in the x-, y- and z-direction, each pair 
yielding one "up" and one "down". This gives us the likelihood function (cf. 0, eq. (9)] 29, 
eq. (4.2)]), 

(which we convert to the spherical coordinates (jSJ) in which we perform our Mathematica 
computations). 

Then, we have S KL (post B \\p Bq=ltrunc ) = 0.169782 and S KL (post BureSq=1 \\p B ) = 0.197657. 
The relative magnitudes of the information gains obtained by passing from priors to posteri- 
ors (0.101846 to 0.169782 and 0.0661775 to 0.197657) seems to suggest that p B is somewhat 
more noniformative than p Bq=l trunc- This is confirmed, using the testing structure given in 
jlll 1^ (cf. jj|), if we formally use a likelihood (L(x,y, z)^), which is the square root of 
(HHJ), to compute post B and post BureSq=1 . Then, we see a decrease in relative entropy from 
0.101846 to 0.093849 and an increase from 0.0661775 to 0.114669. So, p B can be made closer 
to p Bq=1 trunc by adding information to it, but not vice versa, leading us to conclude that p B 
is more noninformative than p Bq=1 trunc, since it assumes less about the data. (Let us note, 
however, that in the class of monotone metrics llgl, the Bures or minimal monotone metric 

n , 

appears to be the least noninformative (cf. |13l . sec. 5]). The maximal monotone metric, on 
the other hand, is not normalizable to a proper prior probability distribution over the Bloch 

is an interesting question of whether there exists a single, distinguished 
normalizable monotone metric which is maximally noninformative.) 



IV. FISHER INFORMATION METRIC OF HUSIMI DISTRIBUTIONS 



Let us now move to a classical context, employing the (generalized) Husimi distributions 
jjy], rather than density matrices to represent the two-level quantum systems. Use of the 
Fisher information (monotone) metric |32j, |33] is now indicated. To generate the (properly 
normalized) escort Husimi distributions {Hr q \) (cf. j^J), from the Husimi distribution (H = 
.£/{i}), we employ the formula (cf. ©), 

-1 



H 



2(r 



qr) 



r) 1+q + (1 + r) 1+q 



H q . 



(17) 
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FIG. 6: Two-dimensional marginal of the four-dimensional extended Husimi volume element (|21|) 

The tangential components of the Fisher information metric for the escort Husimi distri- 
butions (Hr q \) are of the form 



((l + r)f F (t))-\ where |30|, eq. (29)] 



f (t)= (-l + g)(-l+t) 2 (-l+t^) 

qK ' q (1 +t) (1 — q + t + qt — t q — qti — t 1+c > + qt l+q ) ' 1 ; 

In |30l sec. V.D], we succeeded in finding similarly general (for all q) formulas for the 
denominators, but not the numerators, of the radial components. 

In Fig. El we show (having to resort to some numerical integrations, since we lack explicit 
[g-general] expressions for certain of the metric elements) the counterpart to Fig. El for 
the four-dimensional extended Husimi metric. Continuing with our numerical methods, we 
obtain the interesting unimodal curve (Fig. EJ) — the peak being near q = 3.59782, with a 
value there of 0.448488. This portrays the one-dimensional marginal Husimi volume element 
over q (cf. Fig. HJ. In Fig. |S] we show the (quite dimcult-to-compute) one-dimensional 
marginal over r (cf. Fig. |3J). (It appears the upturn near r = 1 may be simply a numerical 
artifact. The difficulty consists in that, in some sense, we have to repeatedly perform 
numerical integrations using results of other numerical integrations. It would be of interest 
to see how the curve changes as the range of q G [|, 500] is modified.) 
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FIG. 7: One-dimensional marginal over q of the four-dimensional extended Husimi volume element 
(|21|). There is a peak near q = 3.59782 



elem 




FIG. 8: One-dimensional marginal over r of the four-dimensional extended Husimi volume element 
l|21|). The upturn near r = 1 may be due to (hard-to-avoid) numerical inaccuracy. 

A. Three-dimensional metric 



For the case q = 1, the (unextended) three-dimensional Fisher information metric over 
the family of Husimi distributions takes the form eq. (2)] 



ds F (p,p + dp) 



2r 3 



Here, 



fp{t) 



L±dr 2 + l(l + r)f F 



1 -rw- 1 



1 + r 



dn . 



t 2 -2tloet- 1 



(19) 



(20) 



which is the limiting case (q — > 1) of (jll)) . To normalize the volume element of this metric 
(Unj) to a prior probability distribution (pf), we divide it by 1.39350989 01 . 
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B. Four-dimensional metric 



In the extended (four-dimensional) case (cf. (|14|)). after having set q = 1, we have, 

^(p,p + dp) 2 = (- - * 16 ; 2 )dq 2 (21) 

2r- (-l + r 2 )logJU , , , , , 2 

H r-= dqdr + ds F (p, p + dp) . 

It 1 

(So, the metric tensor here, in the same manner as in the untruncated extended Bures case 
()1U|). is not fully diagonal. We do not truncate the g-extended Fisher information metric (|21|) 
in any of our analyses.) To normalize its (non-null) volume element to a prior probability 
distribution (p Fq=1 ) over the Bloch sphere, we must divide by 0.24559293. 

V. COMPARATIVE NONINFORMATIVITY ANALYSIS 

We have that S KL {p F \ \p Fq=x ) = 0.229666 and S KL (p Fq=1 1 \p F ) = 0.170145. Further, using 
the likelihood (|16jl. based on six hypothetical measurements to generate posteriors, we ob- 
tain SKL(post F ,p Fq=1 ) = 0.70766 and SKh(post Fq=1 \ \p F ) = 0.0641738. So, the comparative 
noninformativity test, which was initially developed by Clarke leads us to a firm conclu- 
sion that the four-dimensional-based probability distribution p Fq ^ is more noninformative 
in nature than the three-dimensional-based p F . 

Additionally, S KL {p B \ \p Fq=1 ) = 0.148269 and S KL {p Fq=1 \\p B ) = 0.0989669. These are 
converted, respectively, to 0.283218 and 0.0842879 if we replace the first arguments of the 
two relative entropy functionals by posterior distributions based on the (formal) square root 
(L(x, y, z)a) of the likelihood function (|16j). Thus, we can conclude that p Fq=1 is also more 
noninformative than p B - 

Further, S KL (p Bq=ltrunc \\p Fq=1 ) = 0.105463 and S KL (p Fq=1 \\p Bq=ltrunc ) = 0.0914175. 
Again, using the formal square root (L(x, y, z)*) of the likelihood, we obtain changes, re- 
spectively, to 0.245602 and 0.0408236. So, our conclusion here is that p Fq=1 is also more 
noninformative than p Bq=1 t rU nc- We already know from j^] that p B is considerably more 
noninformative than p F . 

Continuing along these lines, S KL {p Bq=1 trunc\\PF) = 0.0191948 and S K L(pF\\pB q=1 trunc) = 
0.0234599 (so the two distributions are relatively close to one another). Using (L(x, y,z)*) to 
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generate posterior distributions, the first statistic is altered (slightly decreased) to 0.0143147, 
while the second statistic jumps to 0.1047772. 

So, assembling these several relative entropy statistics, we have the previously indicated 
ordering of the four priors (The conclusions of the comparative noninformativity test 
appear to be transitive in nature, although I can cite no explicit theorem to that effect.) 



A. Relation to Srednicki's Criterion for Priors 

In Fig. H we show the one-dimensional marginal probabilities of the four prior probabili- 
ties over the radial coordinate r in the near-to-pure-state range r e [.995, 1]. The dominance 
ordering in this plot fully complies with that (0) found by the information-theoretic-based 
comparative noninformativity test. (We note that this ordering is not simply reversed near 
to the fully mixed state [r = 0].) Conjecturally, this could be seen as a specific case of some 
(yet unproven) theorem — perhaps utilizing the convexity and decreasing-under-positive- 
mappings properties 0, p. 35] of the relative entropy functional. 

So, the information-theoretic (comparative-noninformativity) test appears to incorporate 
Srednicki's criterion of "biasedness to pure states" [9J]. (Of course, it would be interesting to 
test the consistency between the comparative noninformativity test and Srednicki's criterion 
with a larger number of priors, as well as in higher- dimensional quantum settings (cf. 3o|).) 
Srednicki does not explicitly observe that increasing biasedness to pure states corresponds 
to increasing noninformativity. He asserts that "we must decide how biased we are towards 
pure states". 

Srednicki focused on two possible priors. One was the uniform distribution over the 
Bloch sphere (unit ball). In |lll . sec. 2.2], we had concluded that this distribution was less 
noninformative than pg, in full agreement with contemporaneous work of Hall Q|. The 
second prior ("the Feynman measure"), which Srednicki points out is less biased to the 
pure states than the uniform distribution, was discussed in 



Neither of the two priors 
analyzed by Srednicki corresponds to the normalized volume element of a monotone metric 

HQ. 
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VI. g-EXTENDED INFERENCE 



In the setting of the g-parameterized escort density matrices Q, the factor ^—f- in the 
likelihood (|16J1 . giving the probability (in the standard three-dimensional Bloch sphere set- 
ting) of one spin-up and one spin-down being measured in the z-direction, would be replaced 
by 

r 2 (l + ^) 2 -(-l + ^)V 



and similarly for the x- and ^/-directions. (For q = 1, we recover 

It would be interesting to ascertain if the volume elements of the extended four- 
dimensional (truncated) Bures and Husimi metrics ( fTUJ) and (JUJ)) could be integrated 
over the product of the Bloch sphere and q e [|,oo] and normalized to (prior) probabil- 
ity distributions. Then, using likelihoods incorporating the form (|22|). one could conduct 
the comparative noninformativity test in a /onr-dimensional setting, rather than only the 
i/iree-dimensional one employed throughout this study. It turns out, however, that the three- 
fold integral — holding q fixed — of the truncated volume element of over the Bloch 
sphere is given by our formula (|12|). Therefore, the four-fold integral of the one-dimensional 
marginal over the indicated product region with q G [|, oo] must diverge. So, to achieve a 
proper probability distribution one would have to truncate q above a certain value. 

Continuing along these lines, we omitted q above 500 (and below q — ~) and normalized 
the volume element of the (truncated) extended Bures metric to a proper probability distri- 
bution. Then, the information gain with respect to such a prior, using L q (z), is 0.0597923 
nats of information, while a single up or down measurement yields 0.134651 nats, and two 
measurements along the same axis giving the same outcome leads to an information gain of 
0.349601. The analogous three (slightly larger) statistics, working in the unextended frame- 
work (where q does not explicitly enter, and is implicitly understood to equal 1), using ps 
as prior, are, respectively, \ — log 3 ~ 0.0680544, and 



8 p F q ({\, 1, 2}, {§, §}, 1) - 7T (-5 + log(64)) - 6 - 12 K 



0.140186, (23) 



6 7T 

(where p F q denotes a generalized hypergeometric function and K w 0.915965594177 is Cata- 
lan's constant) and || — log 5 ~ 0.357229. (We encountered numerical difficulties using 
Mathematica in attempting to extend these analyses to measurements conducted in more 
than one direction, unless we restricted q to a range no larger than on the order of 10.) 
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One might also consider the possible relevance of g-analogs of the Clarke comparative 
noninformativity test, using g-relative entropy (Kullback-Leibler) divergence I37I l3q. 



VII. ^-EXTENDED BURES METRIC FOR HIGHER-DIMENSIONAL QUAN- 
TUM SCENARIOS 

A. Four- Variable 3x3 Density Matrices 



In |35j, we considered an extension of the 2x2 density matrices (J4J) to the 3x3 form 
(by incorporating an additional parameter v ) 



' v + z x — iy^ 



P= 2 



r 2 = x 2 + y 2 + z 2 <v 2 ] 0<v<l, (24) 



2 - 2v 
y x + iy v — z J 

The Bures metric was found there to take the form 

d Bn=3 (p, P + dp) 2 = - (- r ~ V dv 2 + ^—dvdr + +^^dr 2 + -dn 2 ) . (25) 

4 \(1 — v){r z — v z ) r z — v z v z — r z v / 

(So, the tangential component is independent of r, as with JHJ) (cf. H).) Normalizing the 
volume element of ([23)1. we obtain the prior probability distribution : 35J, eq. (18)] 

Pb = 3r2 sin g i (26) 

Att 2 v \/i — f Vf 2 — r 2 

We have calculated that the (five-dimensional) g-extension of this metric has a tangential 
component of the form 

((-r + vf - (r + vY) 2 

4r 2 ((-r + v) q + (r + v) q ) ((2 - 2v) g + (-r + v) q + (r + v) q )' 1 ' 

but have not yet been able to derive simple forms for the other entries of this metric tensor. 

Numerical tests appear to indicate that the volume element of this g-extended Bures 
metric tensor is (also) identically zero. 



B. Abe-Rajagopal Two-Qubit States 

Since our first two attempts above to extend the Bures metric from an n-dimensional set- 
ting to an (n+ l)-dimensional framework, by embedding the q order parameter, have yielded 
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metrics (one of them being (|10j)) with zero volume elements, we were curious as to whether 
or not we could obtain, in some other quantum context, a nondegenerate g-extension of the 
Bures metric. In this regard, we turned our attention to the paper, "Quantum entanglement 
inferred by the principle of maximum nonadditive entropy" of Abe and Rajagopal 7] (cf. 
2, eq. (14)]). 

Their principal object of study is a 4 x 4 density matrix [7J, eq. (32)], being ostensibly 
parameterized by three variables, the order (nonadditivity) parameter q, the g-expected value 
b q of the Bell-CHSH observable and its dispersion a 2 . (Two of the four eigenvalues of the 
density matrix are always equal.) 

We applied the Hiibner formula (jSJ) for the Bures metric to this family of 4 x 4 density 
matrices, considering q as a freely- varying parameter, along with b q and a 2 . Computing the 
3x3 Bures metric tensor, and then setting q — 1, we obtain the metric 

ds AbeRajq=1 (p, p + dp) 2 = Y^ d< l 2 + ( 28 ) 
log(-2 V2 b q + a 2 ) - \og( 2V2b q + a q 2 ) d ^ db + 



8V2 

2 log(8 - a 2 ) - log(-2 V2b q + a 2 ) - log(2 V2b q + a 2 ) 2 

22 aqaa q + 



a JL (db ) 2 + ^ db da 2 + 3 °1 (da 2 ) 2 

-32b 2 q + 4(a 2 ) 2[db « } + l6b 2 -2(a 2 y b ^ + A(-8 + a 2 )(-8b 2 q + (a 2 ) 2 ) { ' 

Here, we have 

c = -4 log(8 - a 2 f a 2 (-8 + a 2 ) + 2 log(-2 V2 b q + a 2 ) log(2 V2 b q + a 2 ) (8 b 2 - a q 4 ) 

(29) 

-log(-2v / 2 6 9 + a f , 2 ) 2 (8b 2 + a 2 (-16 + a 2 ) - 4 V2 b q (-8 + a 9 2 )) 

-log(2v / 2 6 (? + o- (? 2 ) 2 {8b 2 + a 2 (-16 + a 2 ) + 4 V2 b q (-8 + 0-/)) + 

4 log(8-a g 2 ) (-8 + a 2 ) (log(-2 V2 b q + a 2 ) (-2 V2b q + a 2 ^ + log(2 V2b q + a 2 ) (2 V2 b q + o, 2 )) . 

Numerical computations indicate that the volume element of the metric dsAbeRaj q (p, P + dp) 2 , 
for any value of q, is zero. 

In the unextended (two-parameter) case, the nondegenerate volume element (with q = 1) 



is 



(-8+<x q )(-8fc2+( CT 2)2) 



dV A beRaj q =i = - A db q da . (30) 
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FIG. 9: Numerical integration estimates of the Bures volume of the (separable and nonseparable) 
AR two-qubit states, function of q. This \ 
and 1 and appears to be so for all (positive) q. 



AR two-qubit states, as a function of q. This volume is known to be f rs 0.7853981634 for q = \ 



I. q-Invariance of Bures Volumes of Separable and Separable and Nonseparable AR States 



In |39|, it was asserted that for the cases q — | and 1, the associated separability 
probabilities of the Abe-Rajagopal (AR) states were equal to the "silver mean", that is, 
\[2 — 1 « 0.414214 (cf. |4jJ , |4jJ. We have reconfirmed these two probabilities, while also 
finding that the Bures volume of separable and nonseparable states is, in both these cases, 
equal to j ~ 0.7853981634. This also appears to be the Bures volume for all positive q, as 
indicated by the results obtained by numerical integration presented in Fig. El The integrand 
employed (that is, the Bures volume element) was 



dV, 



AbeRaj 



(8 - a q )« 2 (a q - 2V2b q y {2^2b q + cr q ) 



V- ~qj \~ q -V qj \-\--q i ~q/ 

\ Q 4 K 2 - 2 ( 2 (8 _ a f q + _ 2^26,)" + (2V2b q + a q )" 



-^db q da q . 



(31) 

It also appears (Fig. [TUJ) that the Bures volume of the separable (only) AR-states is equal 
0.325323 for all positive q and, thus, the separability probabilities (obtained 



to 



tt(V2-1) 



by taking the ratios) are all simply \pl -1« 0.414214 (that is, the "silver mean"). (The 
numerical integration employed to generate Fig. ^] is more challenging — due to the neces- 
sary imposition of the Peres separability criterion — than to create Fig. El so we could not 
obtain as many significant digits.) 

These g-invariance results stand in interesting contrast to the emphasis of Abe and Ra- 
jagopal "that for q > 1, indicating the subadditive feature of the Tsallis entropy, the entan- 
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0.511.522.53 q 



FIG. 10: Numerical integration estimates of the Bures volume of the separable A two-qubit states, 
as a function of q. This volume is known to be exactly 7r ^ v ^~ 1 ^ m 0.325323 for q = \ and 1 and 
appears to be so for all (positive) q. 

gled region is small and enlarges as one goes into the superadditive regime where q < 1" [?], 
p. 3464 and Fig. 1]. But, in terms of the Bures metric (and others we will see below) the 
measure of the region does not change with q. 

Using the Hilbert- Schmidt metric rather than the Bures, we find that the volume 
of separable and nonseparable AR two-qubit states is equal to for both q = | and 
1 and the volume of separable states is equal to ^= for both these values of q, so the 
corresponding Hilbert-Schmidt separability probabilities are simply ~. (If we employ either 
the Wigner- Yanase [monotone] metric 43] or the arithmetic average [monotone] metric 4l| , 
then, for q = 1, we obtain exactly the same [volume] results as using the Bures metric, 
and for q = | — using numerical rather than symbolic methods in the Wigner- Yanase case 
- quite clearly the same also.) So, it certainly appears that the g-invariance of the total 
and separable volumes of the AR-states is metric- independent. Canosa and Rossignoli 0, 
p. 4] have noted that for the AR-states, the "final maximum entropy density is actually 
independent of the choice of /" , where / is a smooth concave function. 

The Bures separability probability (as well as that based on the Wigner- Yanase metric) 
of the (one-parameter) "Jaynes state" in which (unlike the AR-states) no constraint 

on the dispersion is present (and q is implicitly equal to 1) , is 2 arcsin(v^-i) _ q.271887. 
(Again, note the presence of the silver mean — and implicitly in the very next formula.) 
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The Hilbert-Schmidt separability probability is 

sinhr 1 (2 - V2) + Root [#1 4 - 148#1 2 + 68&, 3] 



0.343602. (32) 



To convert from the AR two-qubit density matrix for q = 1 to that for q = ~, we merely 
have to perform the transformation 

4 (86? + (a?) 2 ) 86i<7? ] 

1 46? + (a?) 2 -8a? +32' 1 46? + (a? ) 2 - 8a? + 32 [ ' 1 ' 

2 2 2 2 2 2 J 

Presumably, there is a (more complicated, in general) transformation between AR-states for 
any pair of distinct values of q. So, in retrospect, the g-invariance of the (Bures, Hilbert- 
Schmidt, Wigner-Yanase and arithmetic average) metric volumes is not so surprising, since 
we are simply working within one family of two-parameter density matrices, the various 
g-manifestations of which can be obtained by suitable reparameterizations. Similarly, the 
null nature of the g-extended Bures metric for the AR-states can be seen in this light. 

2. Trivariate Jaynes state using generalized Bell-CHSH observables 

It would be interesting to extend and analyze the AR-states based on modifications of the 
Bell-CHSH observable (cf. Q sec. 3] Q). In fact, we pursued such a line of investigation, 
using 

B a = 2V2(\$ + )(<5> + \-a\q-)(V-\), (34) 

as the observable, where for a = 1, we recover the Bell-CHSH observable employed by 
Abe and Rajagopal ?[ eq. (6)] (cf. P, eq. (15)]). (We utilized the Jaynes maximum 
entropy strategy Q], implicitly taking q = 1 — so, most precisely, we are extending the 
model discussed by Rajagopal in [lj to incorporate generalized Bell-CHSH observables or, 
alternatively, the Canosa-Rossignoli model to included the dispersion.) Then, the volume 
element of the Bures metric, considering a as a parameter, in addition to the expectation 
b\ and the dispersion ai, was null. Considering, on the other hand, a to be simply a fixed 



constant, the bivariate Bures volume element was of the (non-null) form |w -^^dbidaf, 
where 

C=(8a- a 2 ) (a 2 ) 2 - 4V2(a - l)h (a 2 - 4a) a 2 + (35) 
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and 



D = 16V2(a - l)ab\ - %b\ ((of + 8) a 2 - 3a 2 a + af) . 

For a = 1, we recover ()30)1 . so the associated Bures separability probability is the silver 
mean. For a = 1, using the HS-metric now, the total volume of states is -^m and that of 
the separable states is so the HS separability probability is simply |. 

For the case a = 2, we obtained a result of 0.35368 for the Bures volume of separable and 
nonseparable states, and 0.2000322 for the Bures volume of only separable states, yielding 
a separability probability of 0.566392. The comparable results for the Hilbert-Schmidt case 
(the volume element — independent of b q and a 2 — being Z2a (i+a) db q da 2 ) were and 
^7^, with a separability probability of |. 

Exact integration, then, gave the HS volume of separable and nonseparable states to 
equal, in general, 2y /2(a+a 2 ) ' anc ^ °^ ^ e se P ara ble states — but only for a > 1 — to 
-, so the Hilbert-Schmidt separability probability for a > 1 is simply equal to 



3a- 1 



be = 

I — j^. For a = ~, the HS separability probability appeared to be || ~ 0.578125. 

Then, using the integration over implicitly defined regions feature new to Mathematica 
5.1, we were able to obtain the HS separable volumes, for all (real) values of a, 



5 



-1 < a < 
a = 
10) a = -\ 



\ 



2 

2 ~ 2 v 



■\ + & <a< 1 



V 



16V2 

^ (2y/2 - V2a) 0<a<- 
a>l 
a < -V3 

-v / 2<Q<-lV-v / 3<a<-V / 2, J 



16a 2 (a+1) 

16a 2 (a 2 -l) 

-v / 2a 4 +5v / 2a 2 -v / 2 
16a 2 



(36) 
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FIG. 11: Separability probabilities ()37j) based on the Hilbert- Schmidt metric, as a function of 
the parameter a, for the new class of trivariate Jaynes states employing generalized Bell, Clauser, 
Home, Shimony, Holt (Bell-CHSH) observables (|33|l 



and, dividing by the total HS volume 



■2V2(a+a 2 ) 



), the HS separability probability results, 



\a(a + 1) 



ia(a + 1} 



| (-5 + v 7 ^) a(a + I) 



-\(a - 2)a(a + 1) 



j_ 

4a 



la 



1 

a-1 



(a+l)(a 4 -5a 2 +l) 



■la 



-1 < a < 
a + = 
2a + 1 = v 7 ^ 

0<a<i(-l + v / 5)vi(-l + v / 5)<a<l (37) 
a > 1 

a + V3 < 

-v 7 ^ < a < -1 V -v 7 ^ < a < -V2 



In Fig. ^2 we plot these rather interesting/intricate results. The separability probabilities 
are zero at the isolated points a = —1, 0. For a — > ioo^a maximum of | is approached. We 



also see that the "golden ratio" (or "golden mean" ) [5|, |6( [or its inverse, depending upon the 
definition], v/ ^~ 1 , enters into delineating the different segments over which the separability 
probabilities take different functional forms. (It would seem plausible, although we have 
not conducted a full, detailed analysis that the points at which the functional forms change, 
correspond to separability constraints passing from inactive to active roles, and vice versa.) 

Fake entanglement is avoided in this three-parameter model (3, p. 128] Q> P- 4]. It did 
not appear feasible to directly expand the trivariate (b q , a 1 , a) scenario just investigated to 
a quadrivariate one. 
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FIG. 12: Separability probabilities (|38|) based on the Hilbert-Schmidt metric, as a function of the 
parameter a, for the (Canosa-Rossignoli) class of bivariate Jaynes states employing generalized 
Bell, Clauser, Home, Shimony, Holt (Bell-CHSH) observables (|34|) . The variance af has been set 
to its value 4(1 + in the single constraint (Horodecki) case 

3. Bivariate Jaynes state using generalized Bell-CHSH observables 

ri 

Let us, however, consider a related bivariate (Canosa-Rossignoli-type 2, p. 126]) model, 
in which we set o\ = 4(l + y), the dispersion in the single constraint case eq. (13)]. Then, 
the HS separability probabilities take the form (the HS volume here being 2 ^ 3 "a 2 ~+4a ~^ 

-a -I \ (1-2^) < a < -1 

v/(a - 2)a - vM" + 1) - 1 3a + 2^ < 1 
^/a(a + 1) - a 2a + 1 > V$ 

-a + 2^a(a + 1) - 1 | < a < |(-1 + v 7 ^) . 

We represent this in Fig. El In the ±oo limit, the separability probability approaches 
|. The HS separability probability is zero in the interval a £ [— 1, |]. The point a = 1 
corresponds to the use of the standard (ungeneralized) Bell-CHSH observable (-Bi), and to 
the one-constraint Horodecki model Q]. We can see that the associated HS separability 
probability is the silver mean. 



(38) 



C. Tsallis-Lloyd-Baranger Two-Qubit States 

Tsallis, Lloyd and Baranger have considered a scenario in which the probabilities of 
being in either one of the four states of the Bell basis is given in the form ^' ~" : ' 1 '"' 1 - — - 



4 ' 4 ' 
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and ( l+x+ y +z } gj] (The feasible points (x,y,z) lie in a certain tetrahedron j^, Fig. 3].) 
They a ls0 e m bed the, three-paranteter (two- q ubit) 4 X 4 densi* matrix fl e q . (12)] into 
an (unnormalized) four-parameter 4x4 density matrix p, eq. (14)] by introducing the 
g-parameter. 

Upon its normalization and application of the Hiibner formula (jSJ), we obtained the 
corresponding Bures metric, the volume element of which, in a numerical investigation, 
appeared to be zero, in this four-parameter extended case. So, we have, to this point, yet 
to find any nondegenerate g-extension of the Bures metric (if one is so possible). 

In the unextended three-parameter case, if we employ new coordinates of the form, 

x = 1 — 4 cos6i,y = 1 — 4 sin 6*i cos 9 2 , z = 1 — 4 sin^i sin 6> 2 cos #3, (39) 

then we have simply 

d T LB(p, P + dp) 2 = del + sin Q x 2 de\ + sin 6* sin 9 2 2 d9l (40) 

that is, the uniform metric on the 3-sphere. 

The Bures volume of the (separable and nonseparable) TLB-states is 2 |-, while the Bures 
volume of just the separable states themselves is (thanks to a challenging computation - 



involving a cylindrical algebraic decomposition [cad] [46] - performed by M. Trott) - g n ' . 
(The separable states comprise the cube x,y,z e [— 1, |].) Thus, the Bures separability 
probability QflQ of the TLB-states is (quite elegantly) ^ w 0.27324. 

For the Hilbert-Schmidt metric, the volume of the TLB (separable and nonseparable) 
states is and that of the separable states, , so the HS separability probability is |. 



VIII. RAJAGOPAL-ABE METRIC 

We also investigated the possible application of the "generalized 'metric" (based on the q- 



Kullback-Leibler entropy) 43, eq. (16)] to the two-level quantum systems (j3J) — and seeing 
how it pertains to the family of quantum monotone metrics Jl6j|. For the case q = 1 (which 
should reduce to the Kullback-Leibler symmetrized divergence |47l eq. (6)], our calculations 
yielded that the diagonal elements of the "metric" take the form 

- f - log ( — ) + log W) dr 2 - -r (3 + cos 26 1 + 2 sin 29 1 ) log Wd6\ (41) 

2 V 2 + r / 4 
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-ir(3 + cos 26 2 + 2 sin 2# 2 ) log VP sin 2 Q x dQ\, 

where, as we recall, W = j^. So, it clearly can not possess the form required of a quantum 
monotone metric (cf. Sec. Ill A). However, when we attempted to implement equation (6) 
of ((43]) > bypassing the g-framework, we obtained for the (presumably same?) metric 

1 dr 2 - — \0gWdn 2 , (42) 



1 — r 2 2r 

which does not appear to correspond to a monotone metric. 
IX. CONCLUDING REMARKS 



Naudts 



22j introduced the concept of a 0-exponential family of density operators 



(for which the obvious example is <f>(u) = u q ). He showed that the ^-exponential family of 
density operators, together with a family of escort density operators, optimizes a generalized 
version of the well-known Cramer-Rao lower bound. He assumes that certain Hamiltonians 
are two-by-two commuting. Therefore, the quantum information manifold {pe)e is abelian, 
which "is clearly too restrictive for a fully quantum-mechanical theory" . He suggests further 
work to remove this restriction. 

Abe regarded the order of the escort distribution q as a parameter |48[. He studied 



the geometric structure of the one-parameter family of escort distributions using the Kull- 
back divergence, and showed that the Fisher metric is given in terms of the generalized bit 
variance, which measures fluctuations of the crowding index of a multifractal. 
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